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We study the role played by a Chern-Simons contribution to the action in the QEDg formulation 
of a two-dimensional Heisenberg model of quantum spin systems with a strictly fixed site occupation 
at finite temperature. We show how this contribution affects the screening of the potential which 
acts between spinous and contributes to the restoration of chiral symmetry in the spinon sector. 
The constant which characterizes the Chern-Simons term can be related to the critical temperature 
Tc above which the dynamical mass goes to zero. 
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I. INTRODUCTION 

Quantum Electrodynamics QED^2+i) is a common 
framework aimed to describe strongly correlated systems 
such as quantum spin systems in 2 space and 1 time di- 
mension, as well as related specific phenomena like high- 
Tc superconductivitjiiiSiii. Indeed, a gauge field formu- 
lation of the antiferromagnetic Heisenberg model in d = 2 
dimensions leads to a QED3 action for spinons, see f. i. 
Ghaemi and Senthil^, Morinari^^ and also^. 

We consider the tt-Aux state approach introduced by 
Affleck and MarstoniiS,. The site-occupation of the sys- 
tem by a single particle is enforced by means of a proce- 
dure which was suggested by Popov and Fedotov^iSii^. 

Here we concentrate on the behaviour of the spinon 
mass at finite temperature which is dynamically gener- 
ated by an U{1) gauge field when the action contains a 
Chern-Simons (CS) term. Solving the Schwinger-Dyson 
equation for the spinon field we show below how the crit- 
ical temperature Tc (above which the dynamical mass 
vanishes) is affected by the CS term. 

The outline of the paper is as follows. The first section 
sketches the main steps of the QED3 formulation of the 
2d Heisenberg model. Section lTlll introduces the CS term 
and justifies its presence. In section IIVI we show and 
comment the chiral symmetry restoration of the spinon 
field in the presence of this term. 



where {fj„,fi,cr} are anticommuting fermion operators 
which create and annihilate spinons with cr = ±1/2. The 
projection onto Fock space is exact when the number 
of fermions per lattice site verifies ^ fjo-fi-cr = 1- 

cr = ±l/2 

This can be enforced by using the Popov and Fcdotov 
procedur o^i^i^ which introduces the imaginary chemical 
potential fx = i7r/2/3 at temperature adding the term 
liN to the expression given by Q 



H 



where N — '^fj crfi,>^ counts the number of fermions in 

i,(T 

the spin system. 

In 2d space the Heisenberg Hamiltonian given by 
Eq.(^ can be written in terms of composite non-local 
operators {T^ij} ("diffusons")^'^^ defined as 



If the coupling strengths are fixed as 
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II. FROM THE HEISENBERG INTERACTION 
TO THE TT-FLUX DIRAC ACTION 

Heisenberg quantum spin Hamiltonians of the type 



H — — - ^ JijSiSj 



(1) 



with antiferromagnet coupling {Jij} < can be mapped 
onto Fock space by means of the transformation = 



4^/.,T and 



where 77 is a lattice vector {01,02} in the Ox and Oy 
directions the Hamiltonian takes the form 



V' V — 



+ ^)-f,N (2) 



where i and j are nearest neighbour sites. 

The number operator products {nirij} in Eq.Q are 
quartic in terms of creation and annihilation operators 
in Fock space. In principle the formal treatment of 
these terms requires the introduction of a Hubbard- 
Stratonovich (HS) transformation. One can however 



FIG. 1: Plaquette (□) on a two dimensional spin lattice. Cx 
and Sy are the unit vectors along the directions Ox and Oy 
starting from site i on the lattice 



show that the presence of this term has no influence on 
the results obtained from the partition function. Indeed 
both {rii] and {niUj} lead to constant contributions un- 
der the exact site-occupation constraint and hence are of 
no importance for the physics described by the Hamilto- 
nian |2Jl. As a consequence we leave them out from the 
beginning. 

Using a HS transformation in order to reduce the first 
term in Eq. |(2Jl from quartic to quadratic order in the 
fermion operators and / the Heisenberg Hamiltonian 
reads 



FIG. 2: The contour representation of the energy spectrum 
+ /I = — 2A^cos^(fc2:) +cos^(fcy) for kx and ky be- 
longing to [— f-, f] and showing the presence of the nodal 
points (if, ±-|) where the energy is equal to zero. 
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Then the total flux through the fundamental plaque- 
tte is such that (/)'"^ = tt which guarantees the SU{2) 
symmetry 

Under these conditions the Hamiltonian goes over 
to the TT-flux mean-field Hamiltonian 
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fields 



where {A^j} are the HS auxiliary fields. At this point no 
approximation has been made and Eq.Q is exact. 
The fields Aij can be chosen as complex quantities 
|A|e*'^'J . This parametrization introduces gauge 
ij which are defined on the square plaquette shown 
in figure^and can be decomposed into a mean field part 
and a fluctuation contribution — 0™^ -|- S4>ij. The 
amplitude |A[ too may contain a mean-field and a fluc- 
tuating contribution. In the following we assume that at 
low energy the essential quantum fluctuation contribu- 
tions are generated by the gauge field 6(j)ij and neglect 
the amplitude fluctuations in the sequel. 

The 
that 



's are fixed on the plaquette in such a way 



where (f>™--^ is taken to be constant. 

In order to implement the SU (2) invariance in 
at the level of the mean-field Hamiltonian (j^J we 
follow2i2ii^iiSii& and introduce the configuration 
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fk.cr 
•> fc + TT.cr 
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with 
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+ Acos(f)7fe^,fc„ 
+iAsin(f) 7fc,,fe„+7r 



-iAsin(f) 7fe,,fc„+7r 
- Acos(f)7fc^,fc^ 



where the 7j:'s are defined by 



2 (cos k^ai 



(PFP) 



cosfcj,a2). The eigenvalues of Ti^^p^'' read 



ff 

(P. 
^MF 



M± 2A^cos2(fca;) +cos2(fcj^). 

We are interested in the low energy behaviour of the 
quantum spin system described by Q in the neighbour- 
hood of the nodal points (fc^ = ±^,fcj, ~ where the 

(PFP) _ ,{PFP) \ 



energy gap ( lo 



( + ),fe,cr 



u)] , - ' I vanishes. FigI21 shows 



the contour plot of the energy spectrum Lu^^^^p and lo- 
cates the nodal points. We linearize the energies in the 
neighbourhood of these points. 

Following2ii& the spin liquid Hamiltonian Q at low 
energy can be described in terms of four-component Dirac 
spinous in the continuum limit. The Dirac action of this 
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spin liquid in (2+1) dimensions including the phase fluc- 
tuations S(l)ij around the 7r-flux mean field phase 0™''^ has 
been derived in& and reads 

J I -ia^[(™'^'^ + (l-A)a^9'')]a, 
ip is the 4-dimensional Dirac spinon field 



through a plaquette. This leads to the Maxwell-Chern- 
Simons (MCS) action in Euclidean space 



Se — 



(fr< - -a 
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2«M 



(□5^" + (1 - A)9^a'^) 



(6) 



The implementation of the CS action 
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^ ^ la,ka \ 
f lb,kcr 
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V f2bka / 



where f ^ - and f, ,- ( f ^ - and f„ r ) are fermion cre- 
ation and annihilation operators near the nodal points 
(fjf) ((~f;f)) of the momentum k and indices a and 
b characterize the rotated operators 



f - = — ( f- + f- - 



The first term in (jSJ originates from the U{1) symmetry 
transformation ip e'^^^i/j which generates a gauge field 
o-i-i = df^O. The constant g in |SJl is the coupling strength 
between the gauge field and the Dirac spinons ijj. The 
first term corresponds to the "Maxwell" term —\f^uf^'^ 
of the gauge field where /^^ — d^a„ — d'^a^, A is 
the parameter of the Faddeev-Popov gauge fixing term 
-A {dt'a^f^, Sf"" the Kronecker ^, □ = 5^ + the 
Laplacian in Euclidean space-time. This form of the ac- 
tion originates from a shift of the imaginary time deriva- 
tion dr dr + fJ- where fi is the chemical potential intro- 
duced above. It leads to a new definition of the Matsub- 
ara frequencies of the fermion fields^ ip which then read 

OJF,n = ^F,n - fJ^/i = ^{n + 1/4). 



III. MAXWELL-CHERN-SIMONS ACTION AT 
FINITE TEMPERATURE 

A. Justification and implementation 

As shown by Marston^^, only gauge configurations of 
the flux states belonging to Z2 symmetry (±7r) are al- 
lowed. Hence the flux through the plaquette is restricted 
to (j) = (jf^^ + 5(j) = {0, ±7r}. In order to remove "forbid- 
den" U{1) gauge configurations of the antiferromagnet 
Heisenberg model (0 7^ ±7r) a CS term should be in- 
cluded in the QED3 action in order to fix the total fiux 



qCS 



dT / 



{i'^e'^'-'a^dpa,) (7) 



introduces a new constant k. We show below that this 
constant can be fixed to a definite value. 

From the above action ©, the equation of motion of 
the gauge field in Minkowskian space 



(8) 



leads to a relation between a magnetic field and the CS 
coefficienliiSiSS. li B = dia2 — 92^1 is chosen to be con- 
stant in such a way that the whole system experiences an 
homogeneous magnetic field the equation of motion (jS)) 
of the gauge field becomes 



(9) 



The gauge field afj, is related to the phase 9 (r) of 
the spinon at site r through the gauge transformation 
/f^o- e^^^^^^ff,iT which keeps the Heisenberg Hamilto- 
nian ^ invariant. From the definition of tp gets 
tpp„ e^^^^^tpfa- It is clear that 6 (r) is the phase at the 
lattice site r and that (r) = df^9 (r). Hence the mag- 
netic field B is then directly related to the flux (p through 
the plaquette shown in figure ^ 



<i,j>eO 



^ in) - e (r,)) = g / dl.a 



g I dCl\j.B = gCl\jB = irk 
Id 



(10) 



where fla is the surface of the plaquette and k an integer. 
Here the flux is fixed to be to be equal to {0, ±7r}i2.. 

Hence k can be fixed by the flux through the plaquette 
using equations ((117)) and ©. Defining p = J2 < V'JV'cr > 

as the density of spinon one can indeed rewrite equation 
® 



5P 
' B 



g^U 

TT k 
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where J\f is the number of spinons on the plaquette and 
k is an integer, see (|1U|) . Recalhng that 



number of particle p 27V 

number of flux quanta g\B\/ {27:) k 



(11) 



is the filHng factor of the Landau levelSi in the Quantum 
Hall Effect {QHE) one finally gets 



2-K 



(12) 



The present analysis may suggest that the application 
of a real magnetic field to the spin system could allow to 
detect the presence of spinons through the Quantum Hall 
EffeciiSi. We leave this point for further investigations. 



B. The photon propagator at finite temperature 

In this subsection we construct the dressed photon 
propagator of QED^ with an MCS term at finite temper- 
ature in order to gain information about the interaction 
between spinons i}} and the implication of the CS term 
on the dynamical mass generation. 

Integrating © over the fermion fields ip, the partition 
function of the spin system Z[i/),a] = / V {ij^.a) e^^'^ 
with the action Se given by © leads to the pure gauge 
partition function 



Z[a]= j V (a) e 



S.ff[a] 



(13) 



where the effective pure gauge field action Sef / [a] comes 
in the form 



Seff [a] 



dr I d r< - -a^ 



(□5^'' + (1 - \)d^'d'') 
In det [7^ {d^ - iga^^)] 

(14) 



One can develop the last term in the effective gauge field 
action Sef / [a] into a series and write 



In det [7^ {df^ - iga^)] 



In det G 



■Tr|!GF7"a,|" (15) 



where Gp^{k — k ) — i-^;^Tj^S{k — k ) is the fermion Green 
function in the Fourier space-time with k = (^ujpn,kj, 
hence Gp = -i^^{2T:fl35 (k - k'^ . The first term on 



the r.h.s. of equation ((TB|l is independent of the gauge 
field {a^}. It can be removed from the series since we 
focus our attention on pure gauge field terms. The first 
term proportional to the gauge field n = 1 in the sum 
vanishes since trj^^ — 0. Keeping only second order 
terms in order to stay with gaussian contributions to the 
fluctuations one gets the pure gauge action 



Si% [a] ^ I dr I d^r 







"'^r + ^l.^l.y (2^-^^y (2^ 



(□(5^" + {1- \)^^'^'') 



xtr 



kl 



ra^ih - k")^^j''a, (-(fci - k")) 
k ^ ' 



(16) 



The second term in equation (|16|) has been worked out 
in''. The whole action can be put into the form 



c.(2) 



9_y 



2/9 



d^q 



v(-9) Y 



xaa(-q)\/\%l^ ' a^q) (17) 



where A^°l^ = ,2^J+^2) [q^S^^ - q^q^ - ne^^pqp] 



jhoton pr 

time. The one-loop vacuum polarization terrrs^ reads 



X ^q2yi is the bare photon propagator in Euclidean space- 



H^„ = HaA^^+HsB^,. 

= (fll{q^,)+^2{q^.)) A^^ + n3{q^)B^^ 

where Afj_i, and B^^ are Lorentz invariant tensors given 
in the Appendix and 



ni(9p) = — / dxy^x{l - 

TT 



— sinh /3qy^ x{l — x) 

D{X, Y) 
cos 2'Kxm 



n3(«M) = -5 / dx\og2D{X,Y) 

T^P Jo 

with D{X,Y) ~ cosh (^|3q^yx{l — x)j -I- sin(27rxm) . Here 

the photon momentum = (^ujB.m — ^f) with /i = 

{0, 1, 2}, m is an integer and a = 2g^ the coupling con- 
stant. 

The flnite-temperature dressed photon propagator in 
Euclidean space verifies the Dyson equation 



- A^"^ +n 



(18) 
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The inversion of equation H18(l leads to the dressed pho- 
ton propagator with the CS term at finite temperature 



/\{q' + IlA) {q' 



19 H 



+ (k<z) 



2 



(19) 



IV. "CHIRAL" SYMMETRY RESTORATION 

The couphng of the gauge field to the spinon field 
generates a mass for this field^. Chiral symmetry in four 
dimensions requires fermions to be massless. In this space 
a mass term m^'ip changes sign under chiral transforma- 
tions generated by means of the Dirac matrix 75. Hence 
fermions must be massless in order to keep the action in- 
variant. In three dimensions no real 75 matrix can be de- 
fined. However embedding the (2-|-l)-dimensional space 
into a four dimensional space two types of "chiral" sym- 
metries can be defined from 73 and 7^i22i2^ where 73 
and 75 are 4x4 matrices 



73 



which induce "chiral" transformations e*^^^^ and e'^^^^. 
In (2-1-1) dimensions the algebra is completed by 









I \ 




; 


, 75 = * f 


-1 j 



T3 

-T3 



Tl 

-Tl 



, 7 



T2 

-T2 



where {Ti,i — 1,2,3} are the Pauli matrices and Dirac 
matrices verify 'jf^'yi, + 71^7^ = 2(5^^ in Euclidean space. 

Appelquist et al.^^'^^ showed that at zero temperature 
the originally massless fermion can acquire a dynamical 
mass when the number N "^1) of fermion flavors is 

lower than the critical value Nc = 32 /tt'^. Later MarisSi 
confirmed this result with Nc — 3.3. Since we consider 
only spin- 1/2 systems, N = 2 and hence N < Nc. 

At zero temperature the dynamical mass term is 



renormalized by the CS term 



■(»/16)^J 



m(n—0) 

and even the critical value Nc is affected as Nc = 
Nc [1 -I- (16K/a)^] as shown by Hong and Park^^. 

Here we concentrate on the impact of the CS term on 
the dynamical mass generation and show that chiral sym- 
metry can be restored at finite temperature. An explana- 
tion of the mechanism behind this symmetry restoration 
will also be given. 



A. Effective potential at finite temperature 

In the present theory mass is generated in two differ- 
ent ways. First, as shown earlier, the massless photon 



induces a mass for the spinon through the coupling of 
the two fields^. Second, the CS coefficient gives a mass 
to the "photon" (gauge field a^), tomcs = i^- This can 
be seen from the pure gauge equation of motion for the 
dual field 7^ = ^e^-'Pf,^ 



The massive photon induces the same effect (dynamical 
mass generation) at zero temperatureiSi. 

We show now how the photon mass k (the CS coeffi- 
cient) affects the effective potential at finite temperature 
between two spinous. 

The static effective potential between spinous with op- 
posite charge g is given by 



V{R) = -g^ / drAoo (r, R) 
Jo 

cPq 



27r J (27r)2 
,2 



AooU" = 0,<fe*«-« 



Jo{qR) is the Besscl function of the first kind and 



Ann = 



(g^ + HB(m^O))+ (^.X(».^o)) 

At large distances q the one-loop vacuum po- 
larization parts become H^(m = 0) 9^-1^ ^^'^ 

II_B(m = 0) — -^\n2 where the integer m is related 

to the photon energy, see above H18|l . Hence the longi- 
tudinal part of the photon propagator Aqo leads to the 
definition of a correlation length 



Aon (0, q) 
where ^„ is given by 



nf3 



■In 2 



1 



a/3 
12Tr 



Integrating over the photon momentum q at large dis- 
tance R the effective potential at finite temperature reads 



V{R.I3) 



a 

'n 



SttR'' 



qJojqR) 
q^+S.^' 



which shows that the stronger k the shorter the correla- 
tion length Hence variations of k affects the correla- 
tion length between spinous. Moreover the variation of 
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the flux through the square plaquette also affects the cor- 
relation length since the CS coefficient is related to the 
flux through equation © . If the flux (j) through the pla- 
quette increases the correlation length also increases, 
the larger k the shorter the interaction between spinous. 



B. Dynamical mass generation 

We show how the CS term affects the "chiral" restoring 
transition temperature of the dynamical mass generation. 
The Schwinger-Dyson equation for the spinon propagator 
at finite temperature reads 



G~\k) = 



(27r)2 



(20) 



where p = {po = uJF,mP), G is the spinon propagator, 
Fi, the spinon - photon vertex which will be approxi- 
mated here by its bare value gj^, and A^^ is the dressed 
photon propagator (|19|l . The second term in H2U|) is the 

fermion self-energy E, (G-i = G'-°^ - S). Performing 
the trace over the 7 matrices in equation (|20|l leads to a 
self-consistent equation for the self-energy 



(21) 



In the low energy and momentum limit S(fc) = 
i(/3, k) ~ S(0). Equation simphfics to 



i2n) 



2 "MMV 



1 



-I- m{P, kY 



(22) 



If the main contribution comes from the longitudinal part 
Aoo(0, —P) of the photon propagator H22|l goes over to 



1 = 3-y 



dPP 



(2^) 



p2+nB(m = 0) 



kP 



p2 + n4(m = 0) 



p2 



(23) 



where 




FIG. 3: Dependence of the dynamical mass rn{T) on the 
temperature. Full line : k/A = 0. Long dashed line : 
k/A = 5.10^. Dashed line : k/A = 7.10^. AH curve are 
obtained with a/A = 10^. 



Ua (m = 0) = Hi (m = 0) + Ha (m = 0) 

OlP f ^ 

— / dx^J x{l~ x) tanh PP^J x{l — x) 

I" Jo 

Hb (m = 0) = (m = 0) 

= [ dx In 2 ( cosh |3P^/x(T~■. 
Trp Jo ^ 



Performing the summation over the modified fermion 
Matsubara frequencies ujp.n - 



■-{n + l/4>i& the self 



consistent equation takes the form 



1 = 



(a/A) f' 




As defined above a = g^N with N — 2 since we have 
implemented the Popov-Fedotov procedure^. Here A is 
the UV cutoff and can be identiffed as the inverse spin 
lattice spacing. Equation 124|l can be solved numerically. 

Figure 131 shows the dependence of the dynamical mass 
on the temperature for different value of k/A. This mass 
which is different from zero for low temperature T van- 
ishes at some temperature Tc which depends on n. As 
k/A increases the chiral symmetry transition temper- 
ature Tc/A decreases following the relation y^l^^Q^ = 

g-a(Q/A)-^ where a (a/A) is a coefficient depending on 
a/ A. The behaviour for a fixed a is shown in Figure 01 



7 



n.oi 




1 ' ' ' ' 1 

2000 4000 6000 8000 10000 



FIG. 4: Chiral symmetry transition temperature depend- 
ing on the CS coefficient Here ^ = 10^. 



One can understand the mechanism of chiral symmetry 
restoration as follows. The photon (gauge field) gives a 
mass to the fermions (spinous) through a dynamical mass 
generation mechanism. When the temperature increases 
this mechanism is lowered by fluctuations, the fermions 
gain in mobility (the dynamical mass m {f3, k) decreases). 
This is similar to the situation in plasmas. When the 
temperature is high enough the charged particles com- 
posing the plasma are considered as free particles. Below 
some temperature these charged particles are screened 
thus their mass is renormalized and gets larger than in 
the high temperature plasma. The CS mass k contributes 
also to the photon mass, the interaction V{R, (3) between 
fermions is weakened as k increases for a fixed length 
R, the correlation length gets weaker and thus the 
screening effect gets weaker. Finally the chiral symmetry 
restoring temperature decreases with increasing k since 
the screening effect is smaller and thus fermions gain mo- 
bility, their mass term is renormalized to a smaller value. 
At zero temperature the dynamical mass decreases as n 

increases like "j^^^gj — e (o/is)^^, thus proving that 
the screening effect is lower as the photon mass k in- 
creases and reduces the dynamical mass of the spinon. 

Going back to the Heisenberg model and looking for 
consequences of this chiral symmetry restoration on the 
energy spectrum of spinous, one can see that when the 
CS coefficient differs from zero the gauge field gets a 
gap due to the mass term k. As k goes to oo the spinon 
energy gap m(/3, k) decreases leaving a gapless spinon 
spectrum for a transition temperature Tc going to zero. 

For fixed a/ A the knowledge of Tc fixes k (Fig. 4) and 
consequently the generated mass m{f3,K) (Fig. 3). 



V. CONCLUSIONS 

The low energy spectrum of the Heisenberg model de- 
scribing a two dimensional antiferromagnet quantum spin 
system with an SU(2) symmetry has been mapped onto 
a (2 -|- l)-dimensional quantum electrodynamics action 
with a, U{1) gauge field symmetry^. In this framework 



we showed that the addition of a Chern-Simons term to 
the Maxwell term in the pure gauge field theory of the 
QEDs at finite temperature affects the spinon energy 
spectrum through a chiral symmetry restoration tran- 
sition. The chiral symmetry transition temperature Tc 
above which the dynamical spinon mass is equal to zero 

follows the relation = e""'"/'^^^^ where a (a/ A) 

is a coefficient depending on a/ A and A the UV cutoff. 

The effective potential between two spinons with oppo- 
site charge g at finite temperature shows that the Chern- 
Simons term controls also the screening of this interac- 
tion through the photon mass which is identified with the 
Chern-Simons coefficient k. The correlation length de- 
creases showing that the screening effect gets weaker as 
K increases. 

The value of k can be controlled by fixing the flux 
through a plaquette going around neighbouring spin lat- 
tice sites. The gap in the spinon spectrum shrinks to zero 
with increasing k for a fixed temperature. 

Hence the Maxwell-Chern-Simons term at finite tem- 
perature which is aimed to fix the correct U{1) gauge con- 
figuration provides an interesting way to control the chi- 
ral symmetry restoration temperature and the effective 
potential between spinons. The present study has been 
done in the framework of a non-compact theory. One 
may ask how it will be infiuenced by the presence of in- 
stantons in a compact description of the gauge field. This 
point related to the confinement/deconfinement problem 
is still under discussionSS'2SiSL2£. 
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VI. APPENDIX 

One may believe that a system at finite temperature 
breaks Lorentz invariance since the frame described by 
the heat bath already selects out a specific Lorentz frame. 
However this is not true and one can formulate the sta- 
tistical mechanics in a Lorentz covariant form^. 

We consider a system in 2 space and 1 time dimension. 
Define the proper 3-velocity u^^ of the heat bath. In the 
rest frame of the heat bath the three velocity has the form 
m'' ~ (1, 0, 0) and the inverse temperature (3 characterizes 
the thermal property of the heat bath. 

Given the 3-velocity vector u'^ one can decompose any 
three vector into parallel and orthogonal components 
with respect to the proper velocity of the heat bath, 
the velocity u^. In particular the parallel and transverse 
components of the three momentum with respect to 
read 
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are orthogonal to q'^ 



r = q^ 



111 



(25) 



(26) 



Similarly one can decompose any vector and tensor into 
components which is parallel and transverse to a given 
momentum vector g'^ 



(27) 
(28) 



„ U-u Uji/ 

q^ 



(29) 
(30) 



Since one considers a spin system at finite temperature 
and "relativistic" covariance should be preserved the po- 
larization function may be put in the general formSS 



n„ 



(31) 



It is easy to define second rank symmetric tensors con- 
structed at finite temperature from q'^, u'^ and d^i, which 



and the Dyson equation (|18|) can now be expressed in a 
covariant form if one uses relation H31|) . 
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